Bootstrapping gravity: a consistent approach to energy-momentum self-coupling 
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It is generally believed that coupling the graviton (a classical Fierz-Pauli massless spin-2 field) to its own 
energy-momentum tensor successfully recreates the dynamics of the Einstein field equations order by order; 
however the validity of this idea has recently been brought into doubt 1 1]. Motivated by this, we present a 
graviton action for which energy-momentum self-coupling is indeed consistent with the Einstein field equations. 
The Hilbert energy-momentum tensor for this graviton is calculated explicitly and shown to supply the correct 
second-order term in the field equations; in contrast, the Fierz-Pauli action fails to supply the correct term. 
A formalism for perturbative expansions of metric-based gravitational theories is then developed, and these 
techniques employed to demonstrate that our graviton action is a starting point for a straightforward energy- 
momentum self-coupling procedure that, order by order, generates the Einstein-Hilbert action (up to a classically 
irrelevant surface term). The perturbative formalism is extended to include matter and a cosmological constant, 
and interactions between perturbations of a free matter field and the gravitational field are studied in a vacuum 
background. Finally, the effect of a non-vacuum background is examined, and the graviton is found to develop 
a non-vanishing "mass-term" in the action. 

PACS numbers; 04.20.Cv 



I. INTRODUCTION 

It is a standard view in particle physics that the non-linearity 
of a field theory, such as those of Yang and Mills, can be 
equated with the notion that the field in question carries the 
charge of the very interaction it mediates. This idea has been 
brought to bear on gravity many times, and various arguments 
1121. ISlH H S S l8|] aim to derive general relativity from a lin- 
ear starting point by coupling gravity to the energy and mo- 
mentum of all fields, including the gravitational field itself. 
Despite the conventional wisdom that this self-coupling pro- 
cess is already well understood, Padmanabhan has uncovered 
a number of serious problems with the standard arguments 
jT]. Although we postpone an examination of Padmanabhan's 
analysis to appendix |A] it suffices to express here what is, in 
our view, his most pertinent observation: one cannot start with 
linear gravity, the Fierz-Pauli massless spin-2 action iQl 
and generate the higher-order corrections of general relativity 
by coupling the gravitational field to its own Hilbert energy- 
momentum tensor. More succinctly: one cannot derive the 
Einstein equations by bootstrapping gravi tons' to their own 
energy and momentum. 

To clarify the content of this observation, consider a per- 
turbative expansion of the Einstein field equations Gaji = 
kT^P^'^'^ about a Minkowski background: ^£,^3 = rj^^p + h^p- 
Working to second-order in h^p, we obtain 



where the numbers in parenthesis denote the powers of h^p 

the term contains. Because G^^p = is the equation of mo- 
tion for a massless spin-2 field h^^p, the right-hand side of 
([T]i can be interpreted as this field's source. Thus a satisfying 
physical picture suggests itself: the gravitational field h^p is 
induced by the energy-momentum tensor of all fields T^p = 
tf^p, where t^^p is gravity's own energy-momentum 



12'] 

tensor, identified as ~Gjp/K. In actuality, however, this de- 
scription cannot be formulated in a straightforward manner 
Although the Fierz-Pauli action 5fp is typically used to pre- 
scribe the dynamics of a massless spin-2 field, its Hilbert 
energy-momentum tensor^ 
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is not proportional to G^p, and thus cannot be used as the 
source-term for the second-order field equations. As an alter- 
native approach, one could introduce energy-momentum self- 
coupling at the level of the action: because t^^i^p is a function of 

hf^p, adding the self-coupling term t^ph^P to the Lagrangian 
yields a different result from adding t^p directly to the equa- 
tions of motion. Unfortunately, this procedure also fails to 
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' In discussions of this nature, the word graviton is often used as a short- 
hand for the classical massless spin-2 field. We follow this convention to 
cohere with the literature, but stress that this graviton is in no way quantum 
mechanical. What is actually being referred to is a gravitational wave, a 
classical fluctuation in the geometry of spacetime. 



■ Although other definitions of the energy-momentum tensor exist (see Ct 
we must define according to the Hilbert's prescription j2) in order to 
maintain the analogy with 7"^'^"'"- This definition requires that Spp be "co- 
variantized" (represented in arbitrary coordinates using aflat metric 7„p ) 
and a functional derivative taken with respect to the metric. It is important 
to realise that even though is flat, the arbitrary variations Sy^p required 
to construct the functional derivative inevitably explore cun'ed metrics in 
a neighbourhood of . Thus "covariantization" is not really sufficient: 
the action must be generalised to a curved background spacetime. One 
of the key aims of this paper is to generalise Spp to curved spacetime in 
such a way that energy-momentum self-coupling is consistent with general 
relativity. 



2 



generate —G J^/k in the field equations. 

Padmanabhan claims that these realizations bring to light a 
previously neglected object 5"^ (see appendix |A]i which ap- 
pears to codify the self-coupling of the gravitational field. Un- 
fortunately, this object has many undesirable features: it is 
not a tensor under general coordinate transformations, has no 
clear physical interpretation, and fails to reveal any equiva- 
lence between the coupling of gravity to matter, and gravity to 
itself. 

We propose an alternative solution to this apparent incon- 
sistency: the action for the graviton is not the Fierz-Pauli ac- 
tion but is instead ^2 given by dU, possessing a non-minimally 
coupled term that vanishes when the (vacuum) background 
equations are enforced.^ We shall demonstrate that the 
energy-momentum tensor of this action is the correct second- 
order contribution to the equation of motion, and furthermore, 
that this action provides the starting point for a straightfor- 
ward energy-momentum self-coupling procedure that gener- 
ates the Einstein-Hilbert action (modulo surface terms) to ar- 
bitrary order. We conclude the discussion by extending our 
formalism to non-vacuum spacetimes. 

Throughout the article we employ the abstract index nota- 
tion ifToll . with lower-case Roman indices indicating a tensor's 
'slots', and Greek indices serving to enumerate its compo- 
nents in a particular coordinate system. The metric has sig- 
nature (— , +, +, +), K = 8kG/c'^, and the Riemann and Ricci 
tensor are defined with the following conventions: ^"/^^/V^ = 



II. THE GRAVITON ACTION 

Contrary to the standard approach, we represent the gravita- 
tional field as a perturbation h"^ of the inverse physical metric 
g"'' from the background g"^: 

This expression is exact in that we have not neglected terms 
0{h^y, in contrast, the physical metric gab = gab - h"'gcagdb + 
0{h^) . Following this convention, we use the contravariant 
field h"^, rather than hat, as the fundamental dynamical vari- 
able of the action.^ In general we will write bars over tensors 



^ More precisely, S2 is the action for the graviton in a background spacetime 
with metric in some small neighbourhood of the solutions of the vacuum 
field equations. We use the term vacuum to signify a region without matter; 
this does not necessarily imply the absence of spacetime curvature. 
Any metric theory of gravity will have an ambiguity as to which variable 
g 6 {g"* , gab . y^SS"* ■•■■} should be identified as the true "gravitational 
field". Such a distinction is of no physical consequence and is largely un- 
necessary for a non-perturbative calculation; however for the present dis- 
cussion we are forced to single out a particular field variable for the ex- 
pansion g = g + h. Our aim is to connect gravity to the particle physics 
notion of a spin-2 field and elucidate a simple energy-momentum self- 
coupling scheme that generates general relativity; to this end we are re- 
quired to pick g e {g"'' ,gab} as it is only for these that h is a genuine 
spin-2 field, i.e. a symmetric tensor (not a tensor density) with (lowest- 



derived solely from the background geometry, and adopt the 
usual notational convenience of raising and lowering indices 
with g^'^ and gab-^ 

We posit that the dynamics, energy and momentum of the 
gravitational field /i"'', propagating in a background spacetime 
with metric gab, are all determined (to lowest-order) by the 
following action: 

S^ltX'] ^^J d'xV^h"\G„bcd+H„bcdW, (4) 
where 

Gabcd = 5 {ga{c8d)b - gabgcd)'^^ ~ ^{cgd){a^b) 

+ Jgab^ic^d) + Jgcd^ia^h) (5) 

is a differential operator representing the linearised Einstein 
tensor (see appendix|B]) and 

Hcbcd = jR{gacgdb + jgabgcd)~Rabgcd- (6) 

While Habcd has no obvious geometric interpretation, we in- 
tend to show that its contribution to the action is necessary for 
the consistency of energy-momentum self-coupling with gen- 
eral relativity. Further motivation for this ansatz is given in 
section HU] 

Naturally, if we are to obtain general relativity without at 
first assuming it, we must begin by considering the graviton in 
aflat background spacetime. Nevertheless, we will see from 
the formalism of section |lll] that (provided we use ^2 to de- 
scribe the graviton) energy-momentum self-coupling gener- 
ates the Einstein-Hilbert action even when the background is 
not flat; g"'' need only satisfy the weaker condition 

G„b^Rab-\gabR = 0. (7) 

While this equation expresses the generality of the analysis 
that is to follow, it should be stressed that no knowledge of (|7| 
will be required to assemble the Einstein-Hilbert action order 
by order: a flat background will serve as a perfectly satis- 



order) infinitesimal gauge transfomiation Sh"^ = IV'^e*). Fortunately, it 
is precisely for g E {g"* , g^b } that the necessary energy-momentum self- 
coupling is its most simple: h''''t„b (see WU . These considerations provide 
no criteria for choosing the metric over its inverse as our expansion vaii- 
able, and while this choice only trivially alters the perturbation theory at 
first-order (/i"* ^ —hub) to second-order (the relevant order for 52, t^b, and 
G^j') the two definitions of the ft-field differ by a term of the form h'"^hl'^. 
Our choice of g = g"* is preferable for this article because it simplifies 
the mathematics of the action and energy-momentum tensor. The reason 
for this is explored in ^III El and stems from the fact that any Lagrangian 
for pure gravity must contain more factors of g"* than gcb in order that 
all the derivatives d„ be contracted; thus an expansion in g = will be 
algebraically simpler. Indeed, this observation still holds when coupling 
gravity to a scalar field or a l-formAn, and thus taking g = g"'' simplifies 
many of the calculations of the non-vacuum case also (see i|IVK 
^ The only exception to this rule is the physical metric and its inverse, for 
which g"* ^ gcdf'g"', but rather f'^gbc = 
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factory starting point.^ No matter which background we use, 
however, it is absolutely crucial that we refrain from inserting 
this particular metric (or even equation (|7]i) into the action, 
thereby reducing ^2 to / d'^x^~gh"''Gahcdh"' ■ This is be- 
cause we will need to be able to perform arbitrary variations 
of g"'', not just those consistent with Rabai = or Rah — 0, to 
construct the energy-momentum tensor for h"'^. That said, it 
will be instructive to temporarily ignore this advice so that we 
may relate ^2 to the Fierz-Pauli action. 



A. The Fierz-Pauli action 

For a flat background, Habcd vanishes, and we can choose 
coordinates {x"} such that g"P = 77"^^ and evaluate ^2 as 
a functional of the components h"^. Integrating by parts 
and discarding surface terms, we find that S2 reduces to 
^ / d'^x^pp, where 

+ dahdiih"P (8) 

is the Fierz-Pauli Lagrangian |[l|] J Modulo surface terms and 
an overall rescaling, J^pp is the unique specially relativistic 
Lagrangian for a symmetric tensor field /i"/^ that is invariant 
under the infinitesimal gauge transformation 5h"l^ — Id^'^e^^ 
(see ^ for proof); hence it is the Lagrangian for the graviton 
(massless spin-2 field) in flat spacetime. 

Starting from we can "covariantize" by making the 
replacements rjap ga^, <?« ^ a and multiplying by ^—g. 
This process obviously generates a unique manifestly covari- 
ant Lagrangian density if g"'' is flat, as in this case the pro- 
cedure is equivalent to representing the same Lagrangian in 
arbitrary coordinates. However, for the purposes of calculat- 
ing the energy-momentum tensor (via arbitrary variations of 
g"'^) it will be necessary to generalize ^pp to arbitrary back- 
grounds, and for a curved metric the covariantization proce- 
dure is ambiguous. To see this, observe that we can transmute 
the third term of (O by twice integrating by parts: 

dxh"Pdahp^ ^ dah"Pd^hp\ (9) 

However this equivalence relies on the commutativity of par- 
tial derivatives, and does not occur for the covariant deriva- 
tives of a curved background; instead, integration by parts 



^ Of course, once the self-coupling procedure is complete, and the Einstein- 
Hilbert action has been assembled starting from the graviton on a flat back- 
ground, we will be in a excellent position to justify (7), as this is precisely 
the field equation (applied to the background) that we will have derived. 
With hindsight, then, we can see there was nothing special about our flat- 
space starting point: we may begin with any one solution to and use 
energy-momentum self-coupling to derive the action (and field equation) 
that defines all the others. 

^ Here and elsewhere we use the customary shorthand h = h"^ = h^^g^i,. 



yields 

(10) 

Thus we are forced to make a seemingly arbitrary choice: do 
we to covariantize ([8]) as written, or should we do so after per- 
forming ^1 These two possibilities determine Lagrangians 
which differ by h'"h''^.Rab + h"''h"'Racdb\ they lead to different 
(first-order) equations of motion if the background is curved,^ 
and determine different energy-momentum tensors even if the 
background is flat.^ This last problem is discussed by Pad- 
manabhan |[l[], and is one of his many non-trivial objections to 
the conventional wisdom that general relativity is the unique 
energy-momentum self-coupled limit of the flat-space mass- 
less spin-2 field. 

A greater problem than this ambiguity, however, is that nei- 
ther choice (nor an admixture) leads to general relativity after 
coupling it to its own energy-momentum. As we shall see 
in section [nil the contribution from h"''Hahcdh"' is necessary 
to achieve this, and it is impossible to use the covariantizing 
ambiguity to produce this tensor because it does not contain 
h"'^h"^ Racdb- Instead, the presence of Habcd represents a rather 
different coupling ambiguity faced when moving from a flat 
background to a curved one. Typically we would invoke the 
Einstein equivalence principal to banish from the action terms 
coupling matter fields and Ricci tensors; we would argue that, 
working in locally inertial coordinates about a point p, the La- 
grangian at p should have the same form as the Lagrangian in 
flat spacetime. This amounts to a minimal coupling proce- 
dure: once we have covariantized a specially relativistic La- 
grangian, the job of coupling the field to the gravity is com- 
plete. However, while this rule may make sense to curve the 
background spacetime of a spin-2 field that is "just another 
matter-field" and has nothing to do with gravitation, it is far 
from clear that the principal should hold for the graviton, for 
which it was only ever a convenient fiction to think of as a 
tensor field propagating over a background geometry. 

In summary, the Fierz-Pauli action is insufficient to deter- 
mine S2 for an arbitrary background geometry; the principal of 
equivalence fails to give a unique solution, and cannot justify 
all the contributions necessary for an energy-momentum self- 
coupling procedure consistent with general relativity. How- 
ever, it was never our aim to construct general relativity from 
^pp, and we do not pretend to be able to derive a curved 
spacetime theory of gravity from purely specially relativistic 
concepts. 52 will serve as our starting point, and the only sig- 
nificance we shall ascribe J^pp is that of a special case. 



The first-order field equation only describes the spacetime perturbations 
of general relativity if the ambiguous term is covariantized to become 
Vch"''Vahif; see TO^ and Appendix|B] 
' Note that all other terms of JJpp are invariant under the operation that gen- 
erated (9) so do not introduce further ambiguity. 
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B. Field equations 

Leaving the Fierz-Pauli action behind, we retrain our at- 
tention on ^2 and begin the process of deriving its advertised 
connection to general relativity. First, we shall calculate the 
associated field equations. As usual, the equations of motion 
are derived from the condition that their solutions be station- 
ary configurations of ^2 with respect to variations in the dy- 



namical field h . As we will have no cause to vary g" in 
the derivation, we can enforce the background equations (|7]i 
immediately and discard Habcd- Next, observe that Gabcd is 
"self-conjugate": for any tensor fields A"'' and B"^ 



(11) 



provided either A"'' or B"^ has compact support. Therefore, 
holding g"'' constant and performing a variation Sh"'^ (a sym- 
metric tensor field with compact support) gives rise to a vari- 
ation in the action 



-g5h"^Gabcdh'"' 



(12) 



As Gabcd is already symmetric in its first two indices, we can 
conclude that the equation of motion is 



1 dS, 



^abcd' 



0. 



(13) 



The centrally important feature of this equation is that 



Gabcdh'^" = G^'^b , the linear approximation to the Einstein ten- 
sor under the inverse metric expansion This is particu- 
larly easy to verify for the special case of a flat background 
in Lorentzian coordinates, but is shown to hold more gener- 
ally for vacuum backgrounds in Appendix |B] Thus ^2 pre- 
scribes the correct first-order equation of motion for the gravi- 
ton. In the next section we show that by adding the energy- 
momentum tensor tab of h"'' (determined by S2) to the right 
hand side of ( fT3] l we successfully generate the Einstein field 
equations correct to second-order. 



C. Energy-momentum tensor 

We will now calculate the energy-momentum tensor of the 
graviton and relate it to the second-order contribution to the 
Einstein field equations. We follow Hilbert's prescription and 



define the energy-momentum tensor as a functional derivative 
of the action with respect to the (background) metric: 



tab — 



-1 5S2 
^Sg"'' 



(14) 



where h"'' (rather than hab or /z"^) is to be held constant when 
taking this derivative, as this is the field we have taken to be 
the fundamental dynamical variable. ' ^ 

As an aside, it is worth contrasting the variational definition 
(fl4] l with Noether's (canonical) energy-momentum tensor: 



(15) 



comprising the four conserved currents associated with the in- 
variance of the Lagrangian under rigid spacetime transla- 
tions. The canonical tensor cannot be used in the present dis- 
cussion for a number of reasons. Firstly, it is not uniquely 
determined by the action for h"^: as it depends directly on 
the Lagrangian, we are free to alter f^n by adding a four- 
divergence to without changing either the dynamics of /j"* 
or 82- Secondly, we require a symmetric tensor to act as the 
source for the first-order field equation ( fTSl ). but the canon- 
ical tensor need not have this property.'^ Lastly, Noether's 
definition does not naturally generalize to curved spacetime 
in such a way that t^^a inherits a covariant conservation law 
lITlil . None of these issues arise with tab, ™d in any case 
our aim has been to connect the coupling between matter and 
gravity found in general relativity with a perturbative coupling 
of gravity to itself; it is the Hilbert energy-momentum tensor 
of matter, not the canonical tensor, that appears in the full 
Einstein field equations as the gravitational source. For these 
reasons we discard the canonical tensor and henceforth refer 
to tab, following Hilbert's prescription ( fT4l l. as the energy- 
momentum tensor of /;"''. 

To begin the calculation of tab, we divide the action into two 
pieces ^2 S2G + S2H- 




gh"''Gabcdh^'', 



gh"'Habcdh"'. 



(16) 
(17) 



It will be convenient to perform the functional derivative (fl4] l 
on these two components separately. Focusing first on ^20, 
we integrate by parts'^ so as to remove the second derivatives 



Of course, the resulting field equation will no longer be a stationary config- 
uration of the action S2. In order that this self-coupled equation of motion 
can be derived from the principle of stationary action it will be necessary 
to introduce a third-order correction to the action S3. Naturally, 53 will 
alter the energy-momentum tensor of ft"'' by a term 0(/i^); however, seem- 
ingly by miracle, this will be precisely the third-order part of the Einstein 
field equations. This process continues indefinitely and is explained sys- 
tematically in illlll For the moment we content ourselves with exploring 
the theory to second-order only. 



that it is the energy-momentum contribution from the second-order action 
S2 only. Here we need not make this distinction. 

■ It is true that the canonical tensor can be made symmetric by adding to 



it an identically conserved "correction" (Ja^^'l^'' 
cancels the antisymmetric part of t^^^ . However, 



a function of /i"* that 
if we allow this sort of 
ad hoc adjustment of the energy-momentum tensor, we only exacerbate the 
problem of non-uniqueness. 

' More precisely, one adds to the integrand a divergence of the form 
da{y^[hVhY) = ^/^Vaih^hf that alters Si only by a function of the 
fields on the boundary (or at infinity) and thus may be neglected for the 
purposes of functional variation. 
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from the integrand 
for which we have introduced the abbreviation 

c d 



c d 
(lb ef ' 



(18) 



K. 



ab ef 



"^ga(egf )b - g"'gab8ef " 25^Jf )i^aS^^ 



ba ef 



K ' — K 



(19) 



cib fe 



ef ab 



An infinitesimal variation in the inverse background metric 
5g"^, vanishing on the boundary of the integral, induces a 
variation in the action 



2k 



dK. 



ab ef 



dgPl 



ef 



where 



C" 



hr' {^b8gcd+^c8gbd~^d8g. 



be) 



-i [25;5l,g,)^ - rgbpg,e) ^r5f^ (20) 



is the connection that arises from the variation of the covari- 



ant derivative: 



V + C. We can move the covari- 



ant derivatives off 5g'"^ in the connection term using integra- 
tion by parts, and arrive at an equation of the form 8S2G = 
J d'^x 5gP''[. . .]pq; the tensor density in square brackets is then 
the functional derivative we seek: 

yf^5gP1 2 

( ^Khef'^ _ 1 _ J. e d\ 

\ dgP^ 2^"' 



ab f(p q) 



(p|/k) 



at ef 



K. 



ab f(p q) 



f 



(21) 



Meanwhile, S2H varies by 



5S2H = 



2k 



d^xyf^5R, 



ab 



-h"''h], (22) 



where we have used the background equation (|2| {after the 
variation) to remove the terms proportional to Rdb', these 
would only be significant if we intended to perform further 
variations in the metric. Now, because 



5Rab = 2V[,C'^„„ 

= ihrgapgqb + \dl5l^gp, ~ S;dbgaq) V,V,5^'^^, 



when we (twice) integrate by parts to alleviate 5^"^ of its 
covariant derivatives, we generate a second-order differential 
operator 



Rpqab = \ga(pgq)b^'^ + \gpq^ {ci^ b) - ^ (agb)(p'^ q) : (23) 



72 , 1 



with the property 



j d\^5RabA-" = j d'x^g5gP%qabA"'' (24) 



for all A"^. Therefore, we can conclude from (|22] | that 



K 8S2H _ 1^ 

/=! 5gPi 2 



b{^g'"{^h'+hedh"')-h"'h) 



(25) 



Finally, we have only to combine equations ( |2TI ) and ( 1251 ). 
expand out all the products and derivatives, and assemble the 
outcome into a formula for f„/, as a function of V^h"^. This is a 
straightforward but arduous calculation, and as such we chose 
to complete it with a computer algebra package. The result is 



^tp, = \gpq {hVa^bh"" + 2h'''Va^hh - 2Kby^h"' - hV^h - ^VahV^h - IV^hab^'h'"' + VehJ'Vbh"' 

+ 2VahVbh"^) + \hV^pV^)h - ^hpqVh + \hy\q + K^//\" - {h^^ya^bhpq + {hpqya^bh"^ 

- K(pV'V^)h\ + \hab^(p^q)h"'' - khVaV(ph^^" + IVahV^'hpq + \Vbhap^'h\ - ^Vaflpq^bh"" 

+ l^phab^cih"'' - ^bh\pV^)hJ' - i Vi/zV(p/z ' + \Vbh%Vah\. (26) 
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It is possible to render this formula rather more manageable by working in a gauge with Vah —0,h = 0: 

+ kKb^ipy^^h-' + ^VbhapV'h^ + IVphabV.h"" - Vbh\^y,)hJ' + \Vbh%Vah''^, (27) 



but we will not need this partially gauge-fixed result for this 
present article.''* 

(2) 

Our task now is to compare tab with G^^ and demonstrate 
that the energy-momentum self-coupling of h"'^ (determined 
by ^2) is consistent with general relativity. Details of the cal- 
culation of g|j^' can be found in Appendix iBl the conclusion 

is 



-.(2) 
'ab 



-KU + 0{h^), 



and thus, to second-order, the vacuum Einstein field equations 
are 



^abcd 



h"' = Kt, 



ab 



(29) 



as advertised. 

As a corollary of ( |29] l. we can confirm Padmanabhan's 
observation that general relativity cannot be derived from 
energy-momentum self-coupling the Fierz-Pauli Lagrangian. 
Only once the contribution from Habai is included will Ein- 
stein's gravity result from an energy-momentum self-coupled 
graviton. This realisation casts doubt on Mannheim's recent 
treatment of gravitational energy-momentum lfT2ll . in which a 
tensor is constructed by applying (O to a covariantized Fierz- 
Pauli Lagrangian, rather than 52. 



III. PERTURBATIVE GRAVITY 

Here we develop the formalism to uncover the root cause of 
the second-order energy-momentum self-coupling ( |29] ), and 
reveal how the process continues to arbitrary order. The vast 
majority of this section applies to any metric theory of pure 
gravity'^ and can be generalized to include interactions with 
matter (see ^IVk Only in section HITEl will we commit to gen- 
eral relativity, fix our action S = Seh, the Einstein Hilbert ac- 
tion, and derive the formula ^ for 52. 

We shall concern ourselves with an expansion of the inverse 
metric g"'' about a non-dynamical background which is 



Gauge transformations are covered in i] III CI we note here only that be- 
cause t^i, is not invariant under the infinitesimal gauge transformation 
Sh"'' = 2V'"e*', only the first formula j26t can be used in all gauges. Al- 
though gauge invariance would be a highly desirable property if we in- 
tended to argue that was a physically meaningful tensor in full general 
relativity, it is an impossible request to make of the tensor we seek, which 

(2) 

should be proportional to the gauge dependent tensor G^^l . 

We require only that the dynamics are determined by an action that is a 

coordinate-independent integral of the metric and its derivatives. 



itself an exact solution of the vacuum field equations: 

^ab ^ fb^Xh"\ 

8S[^ 







ab 



(30) 
(31) 



where A, a dimensionless expansion parameter, is constant 
over spacetime. 

Following (l30b . the action of the exact theory S[g\ becomes 
(28) a A-dependent functional of g"^ and h"^, which can be Taylor 



expanded thusly: 



S[g]^S[g + Xh] 



n=0 



where S„ is the "n partial action" given by 



S„[gM^^{d'lS[8 + ^h])l= 



(32) 



(33) 



The derivative acts on each instance of A/i" in the inte- 
grand of S[g + Xh] by Leibniz's law, removing the factor of A. 
The 'bare' left behind may still be covered by spacetime 
derivatives (9„, but these can be moved onto the remainder of 
the integrand by integration by parts. This operation generates 
the usual functional derivative: 



d^S[g + U] 



5g"''{x) 



S[g + lh]. (34) 



In truth, the left hand side of this equation differs from the 
right by the surface term / d'^xdaJ" created when integrating 
by parts. As this is only a functional of the fields on the bound- 
ary (or asx^ if the integral of S runs over the entire man- 
ifold) it will not contribute to equations of motion or energy- 
momentum tensors, the calculation of which are dependent 
only on variations of the field that vanish on the boundary (or 
have compact support). Hence these surface terms may be 
neglected for our present purposes. 

It follows from the repeated application of (|34] | that 



d",S[g + Xh] = 



d'^xh"'' 



5g 



ab 



S[g + Xh], (35) 



and thus the partial actions ( l33l l are given by 



S„[g,h] = ^ 



d^xh"'' 



5g 



ab 



s[g] 



(36) 



An important consequence of this relation is that, using ^2 
as our starting point, we can generate the entire set of partial 
actions {S„ : « > 3} by calculating 



n! 



d'^xh"'' 



§gab 



S2[8,h], (37) 
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which is possible provided ^2 is known in a neighbourhood 
of whichever particular background (a solution of ( [3T1 )) we 
are interested in. Note that the first two partial actions do not 
contribute to the dynamics of So = S[g] is manifestly in- 
dependent of h"'\ and Si vanishes once the background equa- 
tion dSTT i has been enforced. We conclude, therefore, that ^2 
contains all the information necessary to reconstruct the "dy- 
namical" part of the action 



Sdyn[g,h] = Y,^"S„[g,h], 



(38) 



n=2 



which itself contains all the dynamical information of the full 
action S. This is absolutely key to the calculations of sec- 
tion |II] in which we saw the first consequence of this recon- 
struction process, the recovery of the second-order equation of 
motion from an action that one would expect to encode only 
first-order dynamics. 



where, because ^2 is second-order in h" , Gabai will be a linear 
differential operator dependent only on g"'\^'' The equation of 
motion ( l40l ) now takes the form 



^Gabcdh 



cd _ 
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(43) 



«=3 



where it should be taken as given that terms 0{X^^^) have 
been neglected. This is the A^*-order approximation to the 
equation of motion for h"'' that is consistent with the dynam- 
ics of g"'' prescribed by the action S. The first-order contri- 
bution has been separated from the sum so as to evoke the 
picture of a wave equation XGabcdh'^'' = with a source. In 
the next section we will see that the source term on the right 
of ( |43] | is indeed the energy-momentum tensor of the field h"'', 
neglecting terms (9(A^+^ ). 



A. Field Equations 

In general, we could let A be a free parameter and, on de- 
manding 8S[g\/5g'^'^ = for fixed g"'\ derive a A-dependent 
equation of motion Ex [g,h] = for our dynamical field /z"^. 
Any h"'' that solved this equation would correspond to a met- 
ric g"'' = g"'' + Xh"'' that solved the field equations exactly}^ 
However, if we are interested in approximating small varia- 



tions of the metric (i.e. the limit Xh 



ah 



0) we can choose 



some order to which we want the equation of motion to 
hold: 



SS[g] 



= (9(A^+>). 



This is equivalent to 



5h"'' 



= 0{X^+'), 



(39) 



(40) 



where S^^^ is defined by discarding from 5dyn those terms 
that can be neglected in ( [39] l: 



rA'+l 
'dyn 



[gM^j:^"Sn[8M. 

n=2 



(41) 



We shall adopt this 'W*-order approximation" picture for the 
development of our formalism, as we can always write N — °° 
if we wish to discuss the exact theory. 

For the sake of continuity with the previous section, we in- 
troduce the notation 



5S2[g,h] 



gl^ab 



5S[|]/5r''=0 



-gGabcdh"' , 



(42) 



It is advisable to set A = 1 before attempting to solve Ei [g,h] = 0, as this 
constant can always be absorbed into the magnitude of /i"*. Although this 
refinement was convenient for ijll] here we shall keep A as it provides a 
simple method for tracking the powers of /i"* in expressions and is useful 
as a variable for differentiation. 



B. Energy-momentum tensor 

First we shall demonstrate that the dynamical part of the 
action ( |38] | can be generated from ^2 by a simple energy- 
momentum self-coupling procedure. Observe that, as a con- 
sequence of ( |36] |. we have 



5„[^,/.] = -/d^x/.«''^^«-'[^'''] 



5g 



■ab 



(44) 



Defining the rfi^ partial energy-momentum tensor f"^ by ap- 
plying Hilbert's prescription to the n* partial action. 



\,b ■ 



-1 8S„[gM 



we conclude that 



Sn[8,h] 



.n-1 



(45) 



(46) 



This makes manifest the energy-momentum self-coupling 
procedure that allows us to generate the dynamical part of the 
action ( |38] | to arbitrary order, given only ^2. The n* partial 
action is nothing more than the integral of the contraction of 
h"'' with the energy-momentum tensor of the previous partial 
action (divided by — «). The dynamical part of the action is 
therefore given by 



SayVisM = X'%[g,h] 



«=2 



I n+Un 
I 'ab_ 

n + 1 



(47) 



Note that, for the particular case of general relativity (S — 
5eh), the background equation ^ also sets Sq ~ 0, thus 
■^dyn = Ser (modulo surface terms) and the energy-momentum 



The operator Gabcd defined here coincides with the definition in js) once 
S = Seh has been fixed. This is shown in illll El bv deiiving 52 . 
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self-coupling procedure recovers the entire action of the full 
theory, not just the dynamical part. 

Because of factors of n + 1 dividing each f"^ in (07]), it is 
not the case that in the action h"^ couples directly to its (A^*- 
order) total energy-momentum tensor, given by 



^ab — 



'8 



ab 



LA", 



ab- 



(48) 



Instead, the numerical denominators account for the n + 1 fac- 
tors of h"'' in h"^t'^j^, and ensure that the equations of motion 
do indeed have Jj^ as the source. To prove this, note that for 
any symmetric field /"^ (vanishing on the boundary, or with 
compact support) we have 



^,^^ab8Sn[gM 



= /d^x- 



ab 



1=0 



where a = Xjj. 

'dV^ 



111 5h"'' 
^u—Xda- Thus, 



A,ah ^Sn[gM 



§f^ab 



= j,J^XdPaS[g + Xh + al] 
+ ndl-^daS[g + Xh + al' 

= j;rTv.i^ocdr'S[g + U + al] 
= {daS„-i[g + al,h])^^Q 



l=a=0 



X=a=0 



^4^^ab^Sn-l[g,h] 



5g 



:ab 



Hence we have the following important result: 

5Sn[g,h] 5S„-i[g,h] 



gf^ab 

Or, using definition ( l45T l. 

SSn[g,h] 

§fjab 



5g 



ah 



-'Kb'- 



(49) 



(50) 



(51) 



Therefore the equation of motion ( l43T l takes on the form 



XGabcd^' 



cd 



N+\ 

Kk-'Y^xXb'^ 

n=3 



or, recalling 



XGabcd^' 



cd 



ah- 



(52) 



(53) 



We have derived the relation we sought, demonstrating that 
any metric theory of pure gravity can be formulated as a first- 
order wave equation with its own energy-momentum tensor 
as a source. For every > 1, we can derive the equation of 
motion ( l53T l by applying the variational principle to the action 
S^^^\ the left hand side is the wave equation for the linearised 



theory, and the right hand side is the energy-momentum tensor 
prescribed by the action S'^^^. This energy-momentum ten- 
sor is, to some extent, incomplete: it does not include the 
(9(A^+') contribution from the highest-order partial action 
Sn+\- This contribution could be calculated, if so desired, 
and added by hand to the field equations ( |53]) so that the right 
hand side read K"?"^^^', but this equation would no longer be 
a stationary configuration of the action S^^^ . To remedy this, 

we could introduce a correction to the action X^'^^Sn+i that 
would generate the extra term in the equation of motion; the 
appropriate functional is given by ( |46l ) and couples /;"'' to the 
highest-order partial energy-momentum tensor f^^^ But now 
once again the energy-momentum tensor 7"j^^' is incomplete, 
and we can apply this same line of reasoning anew. So long 
as there is no for which f^'^ vanishes identically, this process 
can continue indefinitely, and as N ^ °o the exact field equa- 
tions are recovered, along with the action S^yn — S — Sq — XSi. 

All that remains is to connect our formalism to the specific 
results of the previous section. For the sake of completeness, 
however, we shall first discuss the gauge symmetries of the 
theory, and deduce the conservation law for rj^^ ' . 



C. Gauge transformations 

Because the action S[g] is a coordinate-system independent 
integral, any diffeomorphism : ^ gives rise to a 

gauge transformation of the theory through the action of 0*, 
the map comprising the puUback of on covector indices and 
the pushforward of 0^' on vector indices: 



sWg\^s{g\. 



(54) 



Taylor expanding both sides about g"^ and applying the back- 
ground equation reveals the gauge invariance of the dynamical 
part of the action: 



where 



lab A * „ah 



Xh'"' = (j) 



-ab 



(55) 



(56) 



In the context of an A^*-order approximation, we must insist 
that (p* = 1 +0{X), otherwise these transformations will map 
the small metric fluctuations Xh"^ onto fluctuations compa- 
rable in magnitude to g"''. We can write a general diffeo- 



morphism of this form as 0* = e where is the Lie 
derivative along a vector field ^" — 0(X). The gauge transfor- 
mations of the theory are hence given by 



^ab ^ h'"h = h"'' + 5h"'', 

5h"' ^ A-'f ^^#^r'+E'^^#^/^"', (57) 



n=l 



11=1 



where we have discarded all terms 0{X'^), as these will only 
contribute terms (9(A^+') to the equation of motion, and 
terms (9(A'^+2) to 5^+'. If we wish we can let ^" = e", an 
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infinitesimal vector field, and derive the infinitesimal gauge 
transformation 



{g"'' + ^h"'') N>2, 



(58) 



Because these gauge transformations (infinitesimal or other- 
wise) are symmetries of 5^^', they map solutions of the equa- 
tion of motion ( |53] ) to other solutions. We can therefore use 
the equation of motion to deduce the transformation law for 



^ab- 



cd 



(59) 



This verifies the earlier remark that the energy-momentum 
tensor is gauge dependent, except in the trivial case N — I, 
for which by definition. It may come as a surprise 

that the energy-momentum tensor does not inherit the gauge 
invariance of the action from which it was derived. It should 
be stressed, however, that S^^^ is not identically gauge in- 
variant: the relation (ISST i is only true when the background 
equation is obeyed. For general g"'', the diffeomorphism in- 
variance of S[g] only furnishes the gauge transformation law 
^■^^yV ~ — A55i, the right-hand side of which has a non- 
vanishing energy-momentum tensor responsible for the varia- 
tion in r^^. Equivalently, the gauge dependence of rj^ can be 
seen to result from the non-commutativity of gauge transfor- 
mations and the functional derivative S/Sg"^ used to define 
rj^ [13]; these operations would only commute if the gauge 
mvariance of 5^+' extended to a neighbourhood of the solu- 
tions of the background equation, rather than being confined 
to the solutions themselves. 



infinitesimal vector field e" with compact support. Then ( l60b 
becomes 



0= d^jc 



ds: 



5g" 



■N+l 



5/1"* 



(61) 



Clearly the second term vanishes (to (9(A^+')) if h"'' solves 
the equation of motion i53[ . and thus 



^N+2^ 



= yd^xy^e'^V"7;;j+' + 0(A 
As this equation holds for any e" it follows that 



(62) 



(63) 



is valid up to and including Because this relation 

holds whenever /i"^ solves its equation of motion, and because 
gauge transformations map solutions to solutions, the conser- 
vation law is gauge invariant. 

It is important to recognize that ( |63] | applies to the (A^ + 
1 )*-order energy-momentum tensor: this is the highest-order 
approximation to the energy-momentum tensor that can be 
constructed from our truncated action S^^\ and is a bet- 
ter approximation than the tensor rj^ which features in the 
equations of motion appropriate to this order. Of course, the 
conservation law for follows from ( |63T l by discarding the 
highest-order term, and ensures the consistency of the equa- 
tion of motion (53[ with the identity ^"Gahcdh"' = 0, which 
holds for all /z"^ once the background equation has been en- 
forced. 



D. Conservation law 

It should be expected that 5^^' [8,h] inherits the diffeomor- 
phism invariance of S[g], and that this symmetry endows the 
energy-momentum tensor with a covariant conservation law 
with respect to the background metric. The derivation pro- 
ceeds in close analogy to the proof of V'T^J^""*^'' = from gen- 
eral relativity. 

We again appeal to the diffeomorphism invariance of the 
action ( |54l ) but this time expand S[g] about g"'^ (a solution of 
the background equation) and about (j)*g"'^ (which will 

also be a solution). The result, 

S^,yl'[r8,rh]=S^,;^'[8M, (60) 

affirms that S^^^ is diffeomorphism invariant. Now let be 
an infinitesimal diffeomorphism: 0* = 1 +J!fe for an arbitrary 



Note that diffeomorphism invariance is equivalent to being independent 
of coordinate system, and is a distinct property from gauge invariance as 
defined in OTCl 



E. Constructing tlie graviton action 

It is now time to close the circle of our discussion and con- 
nect the abstract formalism to our earlier calculation. We shall 
derive here the graviton action ^2, the ansatz of section |II] by 
applying the perturbative formalism to the particular case 

- ^ I d'xy^R = Seh [g] , (64) 

the Einstein-Hilbert action. To proceed, we will use equation 
(|36] | to derive Si, and then ^2, by successive functional deriva- 
tives 5/5g"'^ acting on 5eh[^]- The first derivative generates 

Si[gM^^J d^'x^Gahh"' , (65) 

which of course vanishes for all h"^ when g"^ solves the back- 
ground equation Gab = 0. A second variation in g"^ gives rise 
to 

+ Sf'liK.R - hR,a - gcdGabh"'')] . 
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Replacing 5Rah -> 5g^''R, ■dah in accordance with ( |24] |. we de- 
termine 8Si I dg"'' and assemble 

2 J 5^"' 
1 



2k 



'h"'{h,,,R~hR,,i-gcdGabh"')] 



gh"''{Gahcd + Habcd)h 



cd 



(66) 



In the last line we referred to the definitions (|5]l and and 
made use of the identity 



\f-Scd) 



J abed • 



(67) 



This completes the derivation of the graviton action (01 and 
confirms that it can be used as the starting point of an energy- 
momentum self-coupling procedure (|46] | that generates the 
Einstein field equations and the Einstein-Hilbert action (mod- 
ulo surface terms) to arbitrary order. 

The preceding calculation helps to reveal the advantage 
of using /i"*, a perturbation in the inverse metric, as our 
fundamental degree of freedom. Had we instead taken the 
usual approach, expanding gab = gab + ^ hab and taking 
as fundamental, the perturbative formalism would have un- 
folded identically but for the placement of indices. However, 
the calculation of S2 from S'eh would have differed dramat- 
ically. The Lagrangian of S\ would instead be proportional 
to G"^'\)ab, and because the Ricci tensor is naturally covari- 
ant, the variation of G"'' = Rcdg'^'g'"' ~ jRcdg^'g"'' under 5g"'' 
would have been complicated by the extra two factors of ^"^ 
on the first term, compared to the relevant tensor in our ap- 
proach: Gab ~ Rab ~ jRcdg'^'^gab- This trend continues at 
every order; the i)ab convention leads to a greater prolifer- 
ation of terms in each partial energy-momentum tensor be- 
cause the Lagrangian of 5„ has the form (Vo)^(f)„f,)" so must 
be contracted with a further n + 1 factors of g"'' to render it a 
scalar Each of these metric factors generates a term in the 
partial energy-momentum tensor, and thus act as compound 
interest for the process of energy-momentum self-coupling. In 
comparison, our convention leads to Lagrangians of the form 

a)^ {h"^)" , which only need only n—1 additional factors of 
gab-^^ Clearly the inefficiency of the l)ab approach stems from 
the natural covariance of derivative operators ((5„ or V,,) and 
curvature tensors; the advantages of the contravariant expan- 
sion g"'' = g"^ + h"^ are therefore not peculiar to the Einstein 
Hilbert action, and are expected to be even more distinguished 
in higher derivative theories of gravity. 



" There are of course the instances of f^^dcgde in each V^, but these occur 
equally in either convention. 
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This does not mean that all terms in such a Lagrangian will contain only 
n—1 additional factors of gab', there will often be cases in which g"'' is 
contracted with (V^)^ and thus n+l factors of the metric (and its inverse) 
will be present. These cases only represent a small proportion of all possi- 
ble terms, particularly as n becomes large, and are no worse than the terms 
afforded by the f)^;, convention. 



IV. MATTER 

To avoid over-complicating our discussion, we have so far 
focused exclusively on pure gravity. Here we will go some 
way to remedy this simplification, and generalize the formal- 
ism of the previous section to include the perturbations of mat- 
ter fields, and the effects of non-vacuum backgrounds. 

In the most general case, let the action 5 be a functional 
of g"'' and a generic matter field where A will serve as a 
placeholder for any number of internal or spacetime indices. 
We then expand S about a background [g"^ j"^^) as follows: 



= 'i"^ + X\i/^, 

=> s[g,^'] = f^X"Sn[g,h,^,w], 

,1=0 

where g"^ and ^"^ satisfy the background equations 



dg 



ah 



0, 



5S[g,^] 



0. 



(68) 
(69) 

(70) 



(71) 



As before, each partial action can be calculated from the par- 
tial action at the previous order; with matter included, the ap- 
propriate recurrence relation is 



S„ 



where 



-1 5Sn 



hib - 



fk 



^ 5Sn 



(72) 



(73) 



There are two aspects of this coupling scheme that differ from 
pure gravity. The first is immediately apparent: the h"''tab 
term has been joined by an analogous coupling between mat- 
ter fluctuations y/'* and its "source current" Ja- The second 
difference is hidden within the definitions of tab and j^', be- 
cause the {S„} now represent the partial actions for gravity 
and matter together, h"''tab and \j/^jA are no longer just self- 
couplings, and will in general contain terms coupling /i"^ to 
I//'*. In particular, t'^i^ should now be interpreted as the («*- 
order) energy-momentum tensor due to all the fields: h"'', Xj/'^, 
and the background matter 

Proceeding as before, we can now demand that the dynam- 
ical fields /i"* and Xj/^ solve the field equations of the action 
'^dyn' ~ '^n=2 ^"^"^ and generate approximate solutions of the 
exact field equations (prescribed by S) accurate to 0{X^). In- 
stead of using the definition (|42] | for Gabcdi we write the gen- 
eral form of 52, modulo surface terms, as 



2h"%bAW'' + W'''^ABW"), 



(74) 



once the background equations ( ItTI i have been enforced. In 
the above equation, Gahcd^ hibA^ and Wab are linear operators 
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that depend only on background fields, Gahcd and Wab self- 
conjugate, in the sense given by (fTTI) . and I^bA is conjugate to 



P ■ 



(75) 



for all A"'' or provided one has compact support. These 
definitions lead to equations of motion, accurate to 0{X^), as 
follows: 



where 



AG, 



abed' 



ucd 



XWabW" 



KT^h + XKlabAW^, 



(76) 
(77) 



This assumption will mean that the perturbative expansion of 
S can be described by an energy-momentum coupling proce- 
dure only. To see this explicitly, we expand the action about a 
background (^"'',0): 

S[g + lh,X^lf]^Y^r{S^n[gM+S^n[gXW]), (82) 

where each gravitational partial action 

S,n[gM ^ L(dlS,[g + Xh])^^^ 



d^x/i"'' 



lib 



(83) 



rpN V 1 n," j!^ V Q " 

^ab = 2-,^ 'ab^ •' A = 2-, ^ Ja- 



(78) 



Although this formalism is quite general, it is probably too 
general to be usefully employed. Indeed, the complications 
involved in describing matter as a background field and a dy- 
namical perturbation generally serve to obscure the physical 
interpretation of the mathematics. An interesting example of 



this occurs when one tries to rederive V" T 



by apply- 
ing the argument of section HIlD I The result that now follows 



N+l 
ab 



ab 



(79) 



the physical interpretation of which is far from clear Rather 
than continue with this formulation in its full generality, it will 
therefore be more instructive to examine two special cases. 
First, we set = and consider small matter fields A y/'^ 
interacting with Xh"^. Second, by setting y/'* = we can study 
the effect of a background matter field ^'^ on the propagation 
of the graviton. In principal, one could reach these special 
cases starting from the formalism we have just described, but 
it will be simpler and more illuminating to build them up from 
scratch. 



much as before, and the matter partial actions 



S^^n[g,h,^if] = —{^lS^^[g + XKX^|f])^^Q 



- -T(<3i'(A%[^ + A/w])),.o 



1 



{n-2)\ 
1 



dl-^S^[g + XhM)x=a 



5g 



ab 



n-2 



(84) 



Defining the partial energy momentum tensors for h and i/A'^ 



- 1 5Sg„ m„ _ — 1 55ip„ 



Uib 



respectively, we see that the partial actions are coupled as 



(85) 



S„[g,h] = - Jd'^xy 



'ab I ^ab 



gn-1 f'Vn-l" 



ab 



n n — 2 



(86) 



These partial actions lead to the A^*-order equations of motion 



A. Matter perturbations 

In a region where the matter fields are small enough that 
their effects on spacetime curvature can be described by small 
perturbations Xh"'^ in the inverse metric, we can model the dy- 
namics by taking ~ 0, and describe the matter field using 
X\lf^ alone. As it is often the case for gravitational theories, 
let us suppose that the action S is the sum of a gravitational 
action Sg and a matter action S^^: 



S[g,^]=S.[g]+S^[g,^]. 



(80) 



Moreover, for the sake of simplicity, we take ^'^ to be a/ree 
field: 



S^^[g,X^>]^X'S^[g,^>] ^g"\^'^ 



(81) 



XGabcdh"^ 



XWabW = 



^T!:b-L^"('!b+Z") (87) 

n=2 ^ ' 



N r 

E 

n=2 



^X" 



(«-l). 



(88) 



The first equation confirms that the energy-momentum tensors 
of xj/'^ and h"'' combine as the source for the graviton. The 
second equation describes how the coupling between /i"* and 
acts as a source for Note that even when the matter 
field is not free, because Srfi never contains terms linear in the 
matter fields, labA must be at least linear in so we will 
always have I^bA = when 4"^ = 0. 



12 



B. Non-vacuum background 

For a non-vacuum spacetime, we expect to be able to ap- 
proximate (at least to first-order) the behaviour of a gravita- 
tional perturbation by ignoring the perturbations in the mat- 
ter field that it might induce. Alternatively, we may have in 
mind a particular non-vacuum solution of the field equations 
(^^af! ^jjj wish to find nearby solutions (approximate or 
exact) with precisely the same matter content. For these two 
scenarios, we can set i/'^ = and investigate the effect that the 
background 4''* has on the dynamics of h"^ . 

Considerations of this nature highlight an interesting fea- 
ture of our prior discussion of the graviton action. In sec- 
tion [n] we saw the importance of a contribution to the action 
h"''Habcdh^'' that vanished in the vacuum; the obvious ques- 
tion to ask is whether a similar term exists in the non-vacuum 
case, and whether or not it will vanish on the non-vacuum 
background equations. To answer these questions we will de- 
rive the graviton action for a non-vacuum background, which 
will also include the cosmological constant as a special case. 

Let us restrict our attention to general relativity in the pres- 
ence of a matter field: 



(89) 



S^.[g,^'] = 2 U^x^gJ^^{g">\^>^,da'^). (90) 



The factor of two in the definition of the matter Lagrangian 
compensates for our slightly unusual normalization of 
5eh-^' It should be noted that we have assumed that does 
not depend on derivatives of the metric. This is the case for the 
Lagrangians of all the fields of the standard model except the 
spin- J fermion, which in any case should be coupled to grav- 
ity using the vierbein formalism, e.g. 1 14]; such an approach 
is beyond the scope of this article. The results of this section 
can be generalized to allow ^„ to depend on dcg"^ without 
any great difficulty, but this is an added algebraic complica- 
tion that seems to add little insight to our investigation. 

We proceed by expanding the action about a background 
(^"'', ^'^) just as in dMI ) and ( |69l ), but now, as i/'^ = 0, the cou- 
pling scheme ( l72b reverts to the familiar energy-momentum 
coupling of sectionHni Following precisely the same method 
as section nil El we can compute S2 by two successive func- 
tional derivatives (with respect to g"^) applied to The 
first derivative yields 



S\ = 



where 



-1 5S^[g,^>] 



= -2 



dg' 



ab 



(91) 



(92) 



is the energy-momentum tensor of the background matter. 
The second derivative yields the graviton action: 



S2 




h" Gabcdh'' 



(pab-Kfj) h"''h + 2Kh"''h 



abided 



{R + 2k£'^,) (^^habh"' - 



dg"^dg' 



cd 



(93) 



This is the action we sought: the generalization of equation 
dUi to a non-vacuum background. 

If we are only interested in the linear theory, and have no 
wish to calculate the energy-momentum tensor, then we are 
free to enforce the background equation 



Gab — T^T^ibJ 



(94) 



in the graviton action. In sharp contrast to the vacuum case, 
however, the background equation does not reduce ^2 to 
2^ / di^Xy/—gh"''Gabcdh"K or indeed any other covariantiza- 
tion of the massless spin-2 Fierz-Pauli action. Instead, it ap- 
pears as though the background matter has endowed the gravi- 
ton with mass: 



"Gahcdh"' 



a 



(95) 



where the "mass-term" a is given by 



a 



(h"'Kb - ^l^) +Nabcdh"'lfd, (96) 



with 



M^2K[^^-g"''^-^'' 



dg' 



ah 



N, 



abed ■ 



2k 



dg^b^gcd ■ 



(97) 



We refer to a as a "mass-term" because it is quadratic in /z"^, 
free from derivatives, and has been added to the kinetic term 
h"'^Gabcdh"' in the Lagrangian. However, as we will see for 
the specific case of the cosmological constant, a does not by 
itself determine whether the graviton is actually massive, i.e. 
whether it propagates subluminally; the curvature of the back- 
ground will play an equally important role in the field equa- 
tions. In particular, while it is tempting to identify a mass m 
for the graviton according to = M (at least when Nabcd — 0) 
we will soon see that the background matter often sets M <Q, 
so this idea is essentially untenable. 

To explore these issues, it will be instructive to calculate 
a for a few simple examples. First, consider a scalar field 
background <1> with Lagrangian 



All our actions are twice as large as the usual definition. This normalization 
has no effect on the classical equations of motion, but has allowed us to 
define the energy-momentum tensor without a factor of two, simplifying 
the algebra of ^Mml 



^^^-^g"''da^db^^V{^y, 



the mass-term is 



habh"'-W 



(98) 



(99) 
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To ensure that the scalar field has positive energy density, we 
must insist that V (<I>) > 0; hence M < as previously warned. 
Equation ( |99] | can also be used to find the corresponding mass- 
term for a cosmological constant. In this case the Lagrangian 
is = — A/k", which we can reach from by setting 
(9aO = and V = A/jf. Clearly this gives 

aK = ^(hablf^~\h^\ (100) 

which similarly suffers from M < if the cosmological con- 
stant is positive. 

At this point, the reader may be suspicious that the formu- 
lae for and aA (with M < and Nabcd = 0) signify that 
/i"'' is a tachyon in the presence of a scalar field background 
or a cosmological constant. Indeed, if the background were 
flat and M constant over spacetime, we could derive the field 
equations from ( |95] l. observe that their divergence enforces 
the de Donder gauge condition 

d"ha^ - \d^h = 0, 

and, substituting this back into the equations of motion, con- 
clude that the dynamics of the graviton were described by 

{d^^M)h''^ = 0. 

This argument appears to justify the relation np- — M for the 
graviton's mass, and motivate the conclusion that M < be- 
trays tachyonic behaviour It is important to realise, however, 
that the field equation above is of little relevance to the ac- 
tual physical system we were discussing. In reality, M will 
not be constant, and the presence of background matter will 
inevitably preclude background flatness. To understand how 
this last consideration alters the dynamics of the graviton, we 
shall briefly examine the field equation for h"'^ in the presence 
of a cosmological constant. First, we substitute dlOOb into (|95] l 
and derive the field equation 

Gatcdh"'+^{Kb - \ga,h) = 0. (101) 

In contrast to the naive approach, the covariant divergence of 
this equation vanishes identically, and so cannot be used to 
relate V/^/i and V"/;,,^. In place of this, the gauge invariance 
of the vacuum theory remains intact^^, and the field equation 
may be simplified by setting h — 0, V^h"^ = 0: 

V^h,b-2Rd„hy' = 0. (102) 

Surprisingly, the contribution from aA has been cancelled by 
a term proportional to the background Ricci tensor, resulting 



If we wish to extend our discussion of gauge invariance f illll to include 
background matter in general, we would need to account for the gauge- 
fixing implicit in our starting assumption ifr'^ = 0, which is obviously not 
preserved by a (first-order) infinitesimal diffeomorphism Si//'' = Jif^'i'^. 
However, because A is constant over spacetime, no such difficulty arises 
here, and the transformations = —IV^'e^^ remain a symmetry of the 
equations of motion. 



in a field equation that is identical in form to the first-order 
vacuum field equation ( fT3] l in this gauge. Of course, this does 
not indicate that the cosmological constant has no effect on 
the propagation of h"^, only that these effects are limited to 
the constraints imposed on the background geometry by the 
background equation = A^„i,. For this reason, it does not 
seem particularly natural to interpret IRifabc-h'''^ as endowing 
the graviton with a mass; equation ( 11021 ) can instead be under- 
stood as a (partially gauge-fixed) massless spin-2 field equa- 
tion that has been generalised to cosmological backgrounds. 
Quite aside from this, there is also the technical issue of in- 
terpreting the four-index tensor Rated as a mass: only if this 
tensor can be defined in terms of a single scalar variable (and 
the background metric) could the argument be made that this 
single variable described the graviton's mass. For a non-zero 
cosmological constant, the only background with this prop- 
erty is de Sitter space: R^abe = jisdbgae - gdcgab), thus the 
gauge-fixed field equation ( |102| l becomes 

(v^-'^^h"'' = 0. (103) 

If we were so inclined, we might interpret this as a field equa- 
tion for a graviton with = 2A/3, and note that this rela- 
tion has the correct sign for positive A, unlike the formula 
= — 2A suggested by our preliminary inspection of a^. In 
truth, however, further investigation is needed before we can 
either adopt or discard this interpretation. This is not only 
because (I102l i (of which ( 11031 ) is a special case) can be un- 
derstood as a generalisation of a massless field equation to 
cosmological backgrounds, but also because of the subtleties 
involved in interpreting the wave operator in curved space, 
and issues of whether or not to use a conformal coupling. 
Clearly, more work must be done to ascertain the physical 
ramifications of aA. and the "mass-term" a in general, be- 
fore we can understand the degree to which its effects can be 
thought of as giving mass to the graviton. 

Although massive gravitons and the cosmological constant 
were historically viewed as entirely separate concepts, re- 
cent work has brought to light a number of interesting con- 
nections between the two. Deser and Waldron |15] have 
demonstrated that, in (anti-)de Sitter background spacetimes, 
a massive spin-2 field is stable if and only if > 2A/3, or 
m — 0. While it is intriguing that our de Sitter background 
field equation ( |103t suggests precisely the same special value 
of — 2A/3, Deser and Waldron's analysis differs signifi- 
cantly from our own, so this superficial observation may be 
misleading. In particular, whereas our mass-term arises as 
a direct result of the perturbative expansion, Deser and Wal- 
dron add their mass-term to the action by hand. Thus it is 
far from clear that the massive gravitons of their paper corre- 
spond to the physical system considered above. In contrast, 
Novello and Neves flo] claim to prove that = — 2A/3, 
with the implication that A < 0. This approach considers an 
unusual generalisation of the spin-2 field equation to curved 
backgrounds, making a non-standard choice for the covari- 
antization ambiguous term discussed in section 111 Al Thus, 
while their calculations arguably describe a spin-2 field, this 
does not appear to be a natural way to describe the spin-2 field 
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that results from perturbations of the metric (or its inverse) in 
Einstein's theory. It is our intention to disentangle the connec- 
tions between these two approaches, and our own, in a later 
publication. 

For the sake of completeness, we conclude this section with 
an example of a mass-term that can have M > 0, and Nabcd 7^ 
0. Unlike a\, however, we shall not attempt to derive any 
of the implications for the equations of motion. Consider an 
electromagnetic 1-form background A„, with Lagrangian 

^A^-kF^^-lf^f^FacFbd. (104) 

and note that Fat = 2(9[„A^] is independent of the metric. The 
calculation yields 

^-\kF^ [habh"^ - i/l^) - Klf'^h'^^FacFbd, (105) 

which has the aforementioned properties. 



V. CONCLUSION 

Contrary to the prevailing maxim, coupling the classical 
Fierz-Pauli graviton to its own energy and momentum does 
not recreate general relativity order by order. However, there 
is an alternative action for the graviton Q for which energy- 
momentum self-coupling is consistent with Einstein's theory. 
Using this action, the energy-momentum tensor of the gravi- 
ton (|26T l, added as a source to the graviton's first-order equa- 
tion of motion ( fT3] l. builds a field equation consistent with the 
Einstein equations to second-order. Furthermore, the pertur- 
bative formalism developed in section |lll] reveals that our ac- 
tion provides sufficient information to reconstruct general rel- 
ativity to arbitrary accuracy: a simple recurrence relation ( |46] l 
identifies the energy-momentum tensor at one order as the ap- 
propriate contribution to the action at the next. To any order A^, 
this scheme assembles an action that dictates field equations 
(|53] l in which the graviton's A^*-order energy-momentum ten- 
sor is the source. 

The formal machinery used to understand vacuum perturba- 
tions is easily extended to include matter, although the physi- 
cal interpretation of the most general approach, in which mat- 
ter comprises both a background field and a small perturba- 
tion, is less than transparent. Focusing on matter perturba- 
tions separately from non-vacuum backgrounds serves to clar- 
ify the formalism significantly. In a vacuum background, the 
interactions between the graviton and perturbations of a free 
matter field lead to a field equation dSTT i in which the source 
for the graviton is the sum of gravitational and matter energy- 
momentum. This interaction inevitably induces a source in 
the field equations for matter (ISST l. Alternatively, one may 
neglect matter perturbations and examine the consequences 
of a non-vacuum background. In this case, the dynamics and 
energy-momentum of the graviton are prescribed by the action 
( |93] l, generalizing our previous ansatz. Surprisingly, the back- 
ground matter appears to induce a "mass-term" in the graviton 
action, although it is currently unclear to what extent its inter- 
pretation as a mass is valid at the level of the field equations. 



The mass-terms induced by a scalar field ( |99] ), a cosmological 
constant dlOOI l and electromagnetism dlOSt have been calcu- 
lated. 
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APPENDIX A: PADMANABHAN'S ANALYSIS 

The recent article by Padmanabhan [it] unearths many sig- 
nificant shortcomings of the well known arguments BSfilsIl 
that supposedly derive Einstein's equations by coupling the 
Fierz-Pauli graviton to its own energy-momentum tensor. 
Here we attempt to summarize his observations, and explain 
their relation to this present work. 

In broad terms, Padmanabhan's criticisms fall into three ar- 
eas: 

1 . The Einstein-Hilbert action consists of a bulk term (the 

action) and a surface term. The latter includes a 
piece linear in so there can be no way to construct 
it from a self-coupling procedure that starts with an ac- 
tion that is already quadratic in hfj^f^ ?^ 

2. The starting point, the Fierz-Pauli Lagrangian (|8]l, de- 
scribes a Lorentz invariant field theory, and yet the end 
result, general relativity, is generally covariant. It is 
claimed that this metamorphosis only occurs because 
general covariance has been assumed in the various 
derivations, in which case it is "no big deal to obtain 
Einstein's theory". More generally, the classic boot- 
strapping arguments wield ideas developed in general 
relativity (such as Hilbert's definition of the energy- 
momentum tensor) or use knowledge of the end result 
to achieve their goal. Hence they cannot be regarded as 
a derivation of general relativity /rom^rif principles. 

3. The first-order field equation can only take a symmetric 
tensor as its source; the canonical energy-momentum 
tensor (flST l is not necessarily symmetric, and although 
it can be made to be so, this process is not unique. 
Therefore the energy-momentum self-coupling proce- 
dure is ill-defined. The Hilbert definition is uniquely 
determined by the action, but to use it would violate 



The argument given by Padmanabhan is phrased in terms of non-analyticity 
in a dimensionful coupling constant. This form of the argument depends 
on his particular choice of normalization for h^^ and Seh, but is essentially 
equivalent to the statement given here. 
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criticism 2. Crucially, even if we allow ourselves to use 
Hubert's definition, we still fail to recover the correct 
source term for the second-order field equation. 

It is to this very last crucial point that we have devoted the 
bulk of this paper. We now wish to explain our position with 
regards to the first two criticisms, and also Padmanabhan's 
proposed solution to the third. 

1. Our approach expressly avoids discussing surface terms. 
This has greatly streamlined our formalism, and because such 
terms are completely irrelevant for determining field equa- 
tions or energy-momentum tensors, the only price to pay for 
this simplicity is that we can only claim to reconstruct the 
Einstein-Hilbert action modulo surface terms?^ In this sense, 
Padmanabhan's first criticism still stands, although it is un- 
clear whether it has any great importance. If the action is an 
integral over the whole manifold, and asymptotic conditions 
apply to h"'' such that the surface term at infinity vanishes, 
then of course there is no distinction between the Einstein- 
Hilbert action and the action we have constructed. Even if the 
action is an integral over a manifold with a boundary, so long 
as we consider the action to be a functional over all fields with 
a particular boundary configuration (just as we might think of 
the action of a particle as a functional over all paths with par- 
ticular end-points) the two actions differ only by an irrelevant 
constant. Besides, in situations where contributions from the 
boundary really are important, one does not typically use the 
Einstein-Hilbert action anyway: the Gibbons-Hawking-York 
boundary term ITtI [Tsll must be included to remove the de- 
pendence on second derivatives of the metric. This allows the 
field equations to be derived using a variational principle that 
only demands that the variation in the fields (and not also their 
derivatives) vanish on the boundary. 

Padmanabhan's major concern is that the surface term of 
the Einstein-Hilbert action has some quantum mechanical sig- 
nificance. As the nature of quantum gravity has yet to be un- 
derstood, it remains unclear whether or not this is the case. We 
stress once again that the analysis in this paper is purely clas- 
sical, and that we make no claims as to a quantum mechani- 
cal interpretation. Furthermore, it is not even known whether 
the graviton is a useful theoretical object for describing quan- 
tum gravity. We note again that the Gibbons-Hawking-York 
boundary term is usually included in quantum gravity investi- 
gations for which the boundary is not negligible. 

2. It is our view that Padmanabhan's concerns about gen- 
eral covariance are unjustified: we take the position of Wein- 
berg ifTgll . that "general covariance by itself is empty of phys- 
ical content." Any theory (Lorentz invariant or not) can be 
expressed in arbitrary curvilinear coordinates, so the require- 
ment of general covariance cannot, in and of itself, constrain 
the sort of theory one might construct. Rather, the kinemati- 
cal content of general relativity is encapsulated by the equiva- 
lence principle, that the effect of gravity vanishes locally in an 



inertial coordinate system; thus expressing physical equations 
in coordinate invariant notation is an invaluable tool for de- 
scribing how their dynamics are modified by gravity. It is pos- 
sible that when Padmanabhan refers to 'general covariance' 
he is referring to the equivalence principle also. As the latter 
is tantamount to identifying the gravitational field with a dy- 
namical metric, he would certainly be correct to criticise any 
"derivation" that contained such a step; needless to say, we do 
not appeal to the equivalence principle in our approach. 

General covariance aside, though, Padmanabhan's objec- 
tion to the use of curved-space ideas is a valid one, indicating 
that none of the classic arguments constitute a derivation from 
first principles. Our approach certainly makes use of curved- 
space concepts; however our goals are perhaps not quite so 
bold as the other derivations that Padmanabhan has scruti- 
nized: we do not pretend to derive general relativity purely 
from the ideas of Lorentz-invariant field theory. It should be 
stressed, however, that even if some of the kinematical content 
of general relativity is in some way assumed (curved space- 
time, functional derivatives with respect to the metric, etc.) 
it is still a "big deal" to derive the dynamical content of the 
theory, Einstein's equations. 

3. We have already explained our position with regards 
to the definition of the energy-momentum tensor in section 
III CI the only reason that Hilbert's definition is unpalatable to 
Padmanabhan is that his aim is to start with as little curved- 
space mathematics as he can. However, the failure of the 
Hilbert energy-momentum tensor to give the correct second- 
order term for the Einstein field equations is a more signif- 
icant stumbling-block. We have explained our remedy, the 
use of a different starting action, in the body of this paper. 
Padmanabhan, on the other hand, eschews energy-momentum 
self-coupling and introduces a new object S'^^ that he defines 
with the following algorithm. Start with a Lorentz invariant 
Lagrangian ^ {rjap ^ ^'^ap i ^y^ap) expressed in Lorentzian co- 
ordinates {x"}. Replace every instance of rj^p with the metric 

gap to produce a new Lagrangian ^{gap,hap,dyhap); note 
that this is not the same as expressing ^ in an arbitrary coor- 
dinate system because the partial derivatives da have not been 
upgraded to co variant derivatives Vc,. We can now define 



5"'^=2 



ap 



(Al) 



The subscript reminds us that we must set gap = v\ap ^ftei" 
taking the metric derivative, as we are supposedly working 
in Lorentzian coordinates. Padmanabhan claims to be able to 
reconstruct the action by coupling hap to this new object 

5"^. Unfortunately S"'^ has a number of highly undesirable 
properties, suggesting that it is a rather unnatural object, ill- 
defined in its current form.^^ 



Note that this does not nessesarily mean that we have constructed the 
action, only that the integrand of the action differs from \/—gR by some 
total divergence. 



In private communication, Padmanabhan has indicated that he shares our 
concerns about S"^ and does not believe it to be of any fundamental impor- 
tance; hence we present the case against 5"^ for the sake of completeness 
rather than rebuttal. 



16 



Firstly, as it has been constructed from a Lagrangian rather 
than an action, S"l^ depends directly on surface terms. This 
introduces a very large ambiguity, as 5"^ will depend on 
whether we write the integrand of the action in the form {dh)^, 
as Padmanabhan does, in the form hd^h, or as some arbitrary 
combination of both. Each possibility defines a different 5'"'^ 
and (presumably) leads to a different self-coupled limit for the 
graviton. It seems that the only remedy for this ambiguity is 
to artificially stipulate that Jff contain no second derivatives, 
although we note in passing that even this leaves us free to add 
surface terms of the form d"{(j)Aa) in theories for fields other 
than the graviton. 

The second troubling aspect to S"l^ is the "half- 
covariantizing" algorithm used to construct It should be 
clear that this procedure has only been defined in Lorentzian 
coordinates, thus the matrix S"P does not really constitute 
the components of a tensor, as we have not explained how 
their values change when expressed in another coordinate sys- 
tem. There are essentially two ways to extend the definition 
( lAU to include curvilinear coordinates. The trivial solution 
is to construct the tensor S"'' = S"P (daY'^dp)'' using the vec- 
tors {{da)"}, partial derivatives with respect to the Lorentzian 
coordinates used to calculate S"^ in the first place. This ob- 
viously defines a genuine tensor, so the components S" ^ of 
S"'^ in some curvilinear coordinate system {x" } can be cal- 
culated, and they will be related to S"l^ by the usual transfor- 
mation rules. It should be clear, however, that this solution 
is rather unnatural: suppose we have a Lagrangian expressed 
in a curvilinear coordinate system, then the only way to calcu- 
late the components S" in that system is to first transform to 
Lorentzian coordinates, calculate S"^ according to (lAll) . and 
then transform back to our original coordinate system. Also, 
because this process picks out a special set of coordinates, 
there is also no reason to expect that 5"^ can be written as a 
tensorial function of and V„. The natural way to pro- 

ceed would be to generalize the definition (lAll i in such a way 
that we could calculate ^ working in any coordinate sys- 
tem. It might seem that a viable solution would be to define 
the tensor 



ab 



(A2) 



where ^ — ^{gah^hahS chub) is the fully covariant La- 
grangian, and the subscript indicates that the Christoffel sym- 
bols are to be treated as independent of the metric and 
held constant in the derivative. This expression generalizes 
(lAll l to define a tensor S"'^ in a coordinate invariant fashion; 
because the Christoffel symbols are held constant, no term 
arises from a variation of the covariant derivatives, and 5"'' 



The insistence tliat we be able to calculate the components of this object in 
arbitrary coordinates has nothing to do with curved spacetime or general 
relativity. Rather, this reflects the perfectly reasonable expectation that we 
should be able to express Padmanabhan 's self-coupling procedure m flat- 
space spherical polar coordinates, for example, or any other coordinate 
system we choose. 



will reduce to S"'^ in Lorentzian coordinates. This expres- 
sion gives us some insight into the geometrical meaning of 
Padmanabhan's half-covariantized algorithm; in particular it 
reveals that the derivative d/dg^p used to define S"^ is in 
fact exploring geometries (infinitesimally close to Minkowski 
spacetime) with connections that are not metric compatible. 
It is perhaps unsurprising that this F-constant derivative in- 
troduces a new layer of ambiguity to the procedure, as we 
can now alter by adding terms proportional to = ^cgab 
to the Lagrangian. Although this might seem a rather con- 
trived objection, it is in fact a very common consideration. 
For example, suppose the Lagrangian includes a term of the 
form V^/i"^; should we calculate 5"* by acting with d/dg\f- 
on y a{g"'^hch), or should we first commute the metric past the 
covariant derivative, and act on g"'^'Vahcb instead? Note that 
this issue would have been invisible in Lorentzian coordinates 
because 



d^cgef 



dg. 



lib 



-2F 



,{a c-b) 



(A3) 



which we would have automatically set to zero. It seems the 
only way to avoid this uncertainty in S"'^ is to introduce an- 
other artificial constraint on the Lagrangian: we insist that it 
be written in such a way that no derivatives act on the metric. 
This should be achieved by commuting covariant derivatives 
through the metric, rather than integrating by parts, due to the 
aforementioned issues with surface terms. 

We shall take our analysis of 5"^ no further at this time. It 
is still uncertain whether this object can be generalized, nat- 
urally and uniquely, to form a genuine tensor; without such 
a generalization it is difficult to ascertain what sort of mathe- 
matical object the matrix of functions S"^ is supposed to rep- 
resent. Although we cannot claim to have exhausted all possi- 
bilities, the evidence before us suggests, at the very least, that 
this goal is not easily achieved. 

Aside from these technical issues, we should also empha- 
size that, unlike the energy-momentum tensor, S"^ has no 
apparent physical interpretation beyond its supposed role in 
a graviton self-coupling scheme. Energy-momentum self- 
coupling was justified by analogy with matter-gravity cou- 
pling, and advanced by the notion that the energy-momentum 
of all fields should source gravitation. In contrast, the self- 
coupling scheme involving S"^ only serves to set gravity apart 
from the other fields. Furthermore, our solution displays an 
unusual symmetry between the coupling terms in the action 
and source terms generated in the field equations as a result 
(see ijlllBl i; this symmetry is broken by Padmanabhan's self- 
coupling procedure. 



This is the same operation as the derivative used to acquire the Einstein 
equations from the Palatini action | 20], although here we will have no cause 
to perform the complementary derivative d /dY\^. 
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APPENDIX B: EXPANSION OF Gab 

Here we determine the first two terms of the expansion of 
the Einstein tensor 

G„, = Gi;^+Gf,' + 0(/z3), (Bl) 

induced by a perturbation of the inverse metric about a vac- 
uum background: 

^ab ^ fb^h"^^ (B2) 

Gab - 0. (B3) 

The perturbation in the metric is of course fixed by the rela- 
tionship g"^gbc — 5", 

^ gab = gab - hab + hach' + 0{h^). (B4) 

To begin, introduce a connection E"j^^ between the derivative 
operators and V^: 

E"bc = i8"'i^bgcd + ^cgbd ~ ^dghc). (B5) 
This allow the Ricci tensor to be expressed as 

From ( IB5l l it is clear that 

E'f^ = 0, (B7) 
A'c = -\f''{2V(,A)d~^dhbc). (B8) 

+ ir''(2V(,(/i,),r J - V,{Keh',)). (B9) 

Hence the terms of the expansion Rab = R^'^h '^^^ah + 0{h^) 
can be computed as follows: 

^i!' = 2V[,£'i;) (BIO) 
^1^ - 2(V[XS+<yS). (BID 



Thus, 

Gab - ^ab ~ iSabRcd S 

-5^«fo(-Vc.V,,//^ + v2/i), (B12) 

which confirms that Gabcd, as defined in (|5]l, represents the 
linearised Einstein tensor: 

Gabcdh"' ^dp,. (B13) 

In particular, note that both sides of this equation agree on the 
order of the derivatives in V^.V^^/z^j'^; this is the descendant of 
the covariantization ambiguous term discussed in section HTAl 
To find g''^^ , start with 

^(2) _ (2) 1 ( Di:^)3Cd ,Al),rd\ 
Gab - Kb ^2Sab[Rcd8 +Kd'^ ) 

+ ^hahR[llf'', (B14) 

and substitute equations (IB 101 ) and ( IBllI ). followed by dBSl ) 
and ( IB9b . The bookkeeping for this calculation is characteris- 
tically laborious, but is easily accomplished using a computer 
algebra package; the result is 

Gf,' = - Ktah + ^hGabcdh'', (B15) 

where f„/, is given by ( |26] l. As expounded in section HTBI and 
now confirmed by direct calculation ( IB13l l, the first-order ap- 
proximation to the Einstein field equation is Gabcdh^'^ — 0, so 
Gabcdh"' = 0{h^) must hold true at second-order Clearly it 
follows from this that hGabcdh"' — 0{h^), and hence (|28] | is 
verified. 

The third-order difference between G^^ and — Ktab exists 
because the field equation approximated to second-order in 
( |29] l is actually ^J—gG"^/ \/—g — 0; this is of course entirely 
equivalent to the usual form of the Einstein field equation 

Gab - 0. 



[1] T. Padmanabhan, Int.J.Mod.Phys. D 17, 367 (2008), URL 

arXiv : gr-qc/0409089vl 
[2] S. Deser, Gen. Rel. Grav. 1, 9 (1970). 

[3] R. P. Feynman, F. B. Morinigo, and W. G. Wagner, Feynman 

Lectures on Gravitation (Addison- Wesley, 1995), pp. 74-88. 
[4] S. N. Gupta, Phys. Rev. 96, 1683 (1954). 
[5] R. H. Kraichnan, Phys. Rev. 98, 1118 (1955). 
[6] S. Deser, Classical and Quantum Gravity 4, L99 (1987), URL 

http : //stacks . iop . org/0264-9381/4/L99 
[7] b. G. Boulware and S. Deser, Annals of Physics 89, 193 

(1975), ISSN 0003-4916. 
[8] T. Ortin, Gravity and Strings (Cambridge University Press, 

2004), chap. 3.2.' 
[9] M. Fierz and W. Pauli, Proceedings of the Royal So- 



ciety of London. Series A, Mathematical and Physi- 
cal Sciences 173, 211 ( 1939), ISSN 00804630, URL 
|http : / /www ■ j stor . org/ s table/97457 
[10] R. M. Wald, General Relativity (University of Chicago Press, 
1984), p. 437. 

[11] K. Kuchar, Journal of Mathematical Physics 17, 801 (1976), 
URL http : //link . aip . org/link/? JMP/17/801/1 

[12] P. D. Mannheim, Physical Review D (Particles, Fields, Grav- 
itation, and Cosmology) 74, 024019 (pages 6) (2006), URL 
jhttp : //link . aps . org/abstract/PRD/v74/e024019 

[13] G. Magnano and L. M. Sokolowski, Classi- 
cal and Quantum Gravity 19, 223 (2002), URL 
http : //stacks . iop . org/0264-9381/19/223 

[14] A. Lasenby, C. Doran, and S. Gull, Phil. Trans. R. Soc. Lond. 



18 



A 356, 487 (1998), URL arXiv:gr-qc/0405033vl 
[15] S. Deser and A. Waldron, Physics Letters B 508, 347 (2001), 

ISSN 0370-2693. 
[16] M. Novello and R. R Neves, Classical and 

Quantum Gravity 20, L67 (2003 ), URL 

http: //stacks . iop . org/0264- 938 1/2 07L67 | 
[17] J. W. York, Phys. Rev. Lett. 28, 1082 (1972). 
[18] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 2752 



(1977). 

[19] S. Weinberg, Gravitation and Cosmology (Wiley, 1972), pp. 
91-93. 

[20] M. P. Hobson, G. P. Efstathiou, and A. N. Lasenby, General 
Relativity: An Introduction for Physicists (Cambridge Univer- 
sity Press, 2006), chap. 19.10. 



